Evidence for coupling between the Sagittarius dwarf galaxy and 

the Milky Way warp 



Jeremy Bailin 
Steward Observatory, University of Arizona 
933 North Cherry Ave, Tucson, AZ 85721, USA 
jbailin@as . arizona.edu 

ABSTRACT 

Using recent determinations of the mass and orbit of Sagittarius, I calcu- 
late its orbital angular momentum. From the latest observational data, I also 
calculate the angular momentum of the Milky Way's warp. I find that both 
angular momenta are directed toward I m 270°, 6 = 0°, and have magnitude 
2-8 x 10 12 M & kpc km s _1 , where the range in both cases reflects uncertainty in 
the mass. The coincidence of the angular momenta is suggestive of a coupling 
between these systems. Direct gravitational torque of Sgr on the disk is ruled out 
as the coupling mechanism. Gravitational torque due to a wake in the halo and 
the impulsive deposition of momentum by a passage of Sgr through the disk are 
still both viable mechanisms pending better simulations to test their predictions 
on the observed Sgr-MW system. 

Subject headings: Galaxy: disk — Galaxy: kinematics and dynamics — galaxies: 
individual: Sagittarius dSph — galaxies: interactions — Galaxy: halo 



1. Introduction 

The disk of the Milky Way is warped like an integral sign, rising above the plane on 
one side and falling below the plane on the other. This warp is seen both in maps of neutral 
hydrogen (e.g., Diplas & Savage 1991) and in the stellar distribution (Reed 1996; Drimmel, 
Smart, & Lattanzi 2000; Lopez- Corredoira et al. 2002b). The Sun lies along the line of 
nodes of the warp, where tilted outer rings cross the inner plane (see Figure 1). 

Despite their tendency to disperse when isolated (Hunter & Toomre 1969), warps are 
common in external galaxies (Bosma 1981; Briggs 1990; Christodoulou, Tohline, & Steiman- 
Cameron 1993; Reshetnikov & Combes 1998). This has driven many authors to search for 
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universal mechanisms to excite or maintain warps (see Binney 1992, for a review). Many 
of these proposed mechanisms rely on the dark halo to either stabilize warps as discrete 
bending modes within the halo (Sparke & Casertano 1988; but see also Binney, Jiang, & 
Dutta 1998), or to provide the torque necessary to create the warp (Ostriker & Binney 1989; 
Debattista & Sellwood 1999; Ideta et al. 2000). Other proposed mechanisms include the 
infall of intergalactic gas (Lopez-Corredoira, Betancort-Rijo, & Beckman 2002a), magnetic 
fields (Battaner, Florido, & Sanchez-Saavedra 1990), and interactions with satellite galaxies 
(e.g., Huang & Carlberg 1997). 

Each of these mechanisms can, in particular circumstances, produce realistic-looking 
galactic warps. Although no single mechanism appears universal enough to account for all 
warps, the evolution toward a bending mode (even when no discrete mode exists) appears 
enough like an observed warp (Hofner & Sparke 1994) that warping may be a generic response 
of disks to the individual perturbations they experience. In this case, we should look at 
individual warped galaxies for specific evidence of particular perturbations that explain their 
warps rather than search for a universal mechanism that may not exist. 

The Magellanic Clouds have been proposed as the perturbation responsible for the 
Milky Way's warp. While Hunter & Toomre (1969) found that the tidal distortion from 
the clouds alone is not sufficient to cause the observed warp, Weinberg (1998) proposed 
that orbiting satellites could set up wakes in the Milky Way's halo which could provide the 
necessary torque. Tsuchiya (2002) performed self-consistent simulations of such a system 
and confirmed that for a sufficiently massive halo (2.1 x 10 12 M ), the magnitude of the 
torque can be increased enough to cause a warp of the same magnitude as the Milky Way's. 

The Magellanic Clouds orbit about the center of the Galaxy in a direction orthogonal to 
the line of nodes, i.e., near the line of maximum warp (see Figure 1). Garcia- Ruiz, Kuijken, 
& Dubinski (2000) demonstrated that the warp caused by a satellite will have its line of nodes 
oriented along the satellite's orbit. A simple way of understanding this result is to recognize 
that a torque is a transfer of angular momentum, and therefore the disk will acquire angular 
momentum along the same axis as the orbital angular momentum of the satellite which is 
providing the torque, and tilt toward that axis. Therefore, the Magellanic Clouds are a bad 
candidate for producing the Milky Way warp. 

The orbital plane of the Sagittarius dwarf galaxy (Ibata, Gilmore, & Irwin 1994) does 
intersect the line of nodes, suggesting that it may be a good candidate for producing the 
Milky Way warp (Lin 1996). It is located behind the Galactic bulge and is on a nearly polar 
orbit (Ibata et al. 1997). Ibata & Razoumov (1998) performed simulations which suggest 
that the passage of a sufficiently massive Sgr (5 x 10 9 M ) through the disk could produce 
a warp. Alternatively, its gravitational tides or the tides of a wake it produces in the dark 
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Fig. 1. — Schematic drawing of the plane of the Milky Way, the Galactic center (GC), the 
line of nodes of the warp, and the orbits of Sagittarius (Sgr) and the Large Magellanic Cloud 
(LMC). The plane of Sagittarius's orbit intersects the line of nodes and is orthogonal to the 
plane of the LMC's orbit. Not to scale. 
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halo could exert a warp-inducing torque on the disk. 

If Sgr is responsible for the warp, its angular momentum will be coupled to that of 
the warp. In this Letter, I calculate the orbital angular momentum of Sgr, along with the 
component of the Milky Way disk's angular momentum which does not lie in the common 
plane. I show that they have the same direction and the same magnitude. As there is no a 
priori reason to expect them to be within orders of magnitude of each other, this is evidence 
that Sgr is coupled to the warp, and therefore responsible for it. 

2. The angular momentum of satellite galaxies 

The position, distance, and motion of Sagittarius are given in Table 1, along with 
estimates of its mass and orbital angular momentum. The angular momentum can range 
between 1.7 and 8.6 x 10 12 M kpc km s" 1 and is directed toward I = 276°, b = 0°. 

The major uncertainty in this calculation is the determination of the mass. Ibata & 
Lewis (1998) argued that in order for the satellite to have survived to the present day, it 
must have a massive extended dark matter halo and a total M/L ~ 100 in solar units (Ibata 
et al. 1997). However, Helmi & White (2001) found viable models with more moderate 
masses ranging from 4.66 x 10 s M Q for a purely stellar model to 1.7 x 10 9 M for their 
model with an extended dark matter envelope (see also Jiang & Binney 2000, who find that 
if the original mass of Sgr was large enough for dynamical friction to be important, the 
majority of the mass would have been stripped off after a Hubble time leaving a current 
mass of 1-3 x 10 9 M ). The properties of Helmi & White (2001) 's models seem most in 
agreement with the expected properties of dwarf spheroidal galaxies, and therefore I adopt 
0.4-2.0 x 10 9 M as the range of possible masses of the Sagittarius dwarf. 

Table 2 shows the magnitude and direction of the angular momenta of the Galactic 
satellites with measured proper motions, along with that of the Milky Way warp which 
is calculated in Section 3. The orbital angular momentum of the Large Magellanic Cloud 
(LMC) was calculated using data from Kroupa & Bastian (1997). For the remaining satel- 
lites, the mass was taken from Mateo (1998) and the velocity vector from the tabulated 
reference. 



3. The angular momentum associated with the Milky Way warp 

I calculate the component of the disk angular momentum which is due to the warp in 
the Milky Way's disk, i.e., that which is not directed toward the North Galactic Pole (NGP). 
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Table 1. Properties of the Sagittarius dwarf 

Ibata et al. (1997): 

Galactic coordinates 
Galactocentric distance 
Space motion (U, V, W) 
Galactocentric radial velocity 
Galactocentric tangential velocity 

Derived angular momentum: 

Assumed mass (10 9 M ) 0.4 2.0 

Angular momentum (10 12 M kpc km s _1 ) 1.7 ± 0.6 8.6 ± 3.4 
Direction I = 276°, 6 = 0° 



/ = 5.6°, b = -14° 
16 ± 2 kpc 
(232,0, 194) ±60 km s" 1 
150 ± 60 km s" 1 
270 ± 100 km s" 1 



Table 2. Angular momenta of the Milky Way warp and some Milky Way satellites 



Angular momentum (M kpc km s 1 ) Direction Reference 



Milky Way warp 


1.7-8.6 x 10 12 


I 


= 270°, b = 


0° 




Sgr dSph 


1.6-7.3 x 10 12 


I 


= 276°, b = 


0° 


1,2 


LMC 


2 x 10 14 


I 


= 184°, b = 


9° 


3 


Fornax dSph 


3 x 10 12 


I = 


-- 106°, b = - 


-16° 


4,5 


Ursa Minor dSph 


3 x 10 11 


I 


= 213°, b = 


9° 


4, 6 


Sculptor dSph 


1 x 10 11 


I 


= 226°, b = 


7° 


4,7 



References. - (1) Ibata et al. 1997; (2) Helmi & White 2001; (3) Kroupa & Bastian 
1997; (4) Mateo 1998; (5) Piatek et al. 2002; (6) Schweitzer, Cudworth, & Majewski 1997; 
(7) Schweitzer et al. 1995 
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Table 3. Disk parameters 



Model 
12 3 4 



Dehnen & Binney (1998): 



Stellar disk scale length R^* (kpc) 


2.0 


2.4 


2.8 


3.2 


ISM disk scale length Rd,isM (kpc) 


4.0 


4.8 


5.6 


6.4 


Surface density at solar circle S (M pc~ 


- 2 ) 43.3 


52.1 


52.7 


50.7 


Derived warp angular momenta: (10 12 M 


kpc km s 








Stellar disk 


1.10 


2.16 


3.36 


4.73 


ISM disk 


0.52 


1.10 


1.78 


2.57 


Total 


1.62 


3.26 


5.14 


7.30 


Direction 


I = 


270° ± 


10°, b 


= 0° 
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If the disk rises a height h(R) above the plane at cylindrical radius R, then the total angular 
momentum in the disk which is due to the warp is 

L w = [ 2nR 2 vJl(R) , ffl dR. (1) 

J Rw 1 Vw + ^ 2 

The mass distribution of the disk is taken from Dehnen & Binney (1998). The disk 
surface density for a given component in these models is given by 

E(i2)=E d exp(-^--|), (2) 

where S d is the normalization, R d is the scale length of the component, and R m is introduced 
to allow the ISM to have a central depression 1 . R m = 4 kpc for the gas disk and R m = for 
the stellar disk. The relative contributions to the surface density at the solar circle So are 
0.25 for the ISM and 0.75 for the stars. Dehnen & Binney (1998) distinguish between thin 
and thick disk components of the stellar disk, but because these only differ in vertical scale 
height, which does not affect the angular momentum, I treat them as a single component. 
Their models 1-4, which differ primarily in disk scale length, R d) are all acceptable fits to the 
observations, and therefore provide a reasonable range of mass distributions with which to 
estimate the angular momentum. Table 3 gives the essential parameters for the four models. 

The circular velocity, v c , of the disk from 3 kpc to the solar circle is ~ 200 km s _1 (e.g., 
Merrifield 1992). While most measurements at R > R show a rising rotation curve, Binney 
& Dehnen (1997) argue that a constant rotation curve is consistent with the data when the 
correlations between errors are taken into account. I adopt v c = 200 km s _1 at all radii. The 
uncertainty in the angular momentum due to uncertainties in the mass models dominates 
over any error in the circular velocity. 

The height of the warp above the plane as a function of radius, h(R), appears to differ 
for the stars and for the gas. Drimmel et al. (2000) fit Hipparcos measurements of OB stars 
and find 

h ( R) = i( R - R ™) 2 / R h R>R W) (3) 

t0 R < R w 

with the warp starting at R w = 6.5 kpc and scaled by R h = 15 kpc. Binney & Merrifield 
(1998) approximate the m — 1 mode of the ISM warp as 

h(R) = (^ R ~ Rw)/a R>R W (4) 

y 1 10 R < R m v ; 



1 Note that equation (1) of Dehnen & Binney (1998) has a typo which is fixed above (W. Dehnen 2002, 
private communication) 
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where R w = 10.4 kpc and a = 5.6 when converted to i? = 8 kpc as assumed in the Dehnen & 
Binney (1998) models (Tsuchiya 2002). Binney & Merrifield (1998) also fit an m = 2 mode, 
but the net angular momentum of any even m mode is aligned with the angular momentum 
of the flat disk, so it will not contribute. 

I use equation (3) for the stellar disk and equation (4) for the gas disk. The results 
are shown in Table 3. The majority of the angular momentum is contained in the range 
10 ~ R ^ 25 kpc in all models. The Sun lies within 10° of the line of nodes, so L w is 
directed toward 260° < I < 280°, 6 = 0°. 

4. Discussion 

Table 2 shows the magnitude and direction of the angular momentum of the Milky Way 
warp and of the Galactic satellites with measured proper motions. Both the magnitude 
and direction of the angular momentum of Sagittarius are strikingly similar to that of the 
Milky Way warp. There is no a priori reason to expect this; the angular momenta of the 
other satellites with known orbits span three orders of magnitude and almost 180° of galactic 
longitude (although there is a strong tendency for the satellites to have polar orbits with low 
values of \b\, as suggested by Lynden-Bell (1976) and noted in the anisotropic distribution 
about the Milky Way by Hartwick (2000)). The coincidence of the two angular momentum 
vectors is probable evidence that they are dynamically coupled, i.e., that Sagittarius is the 
perturber responsible for the Galactic warp. 

There are three possibilities for the nature of the coupling. The first is a direct gravi- 
tational tidal torque by the satellite itself (Hunter & Toomre 1969), the second is the gravi- 
tational torque of a wake in the Galactic dark matter halo (Weinberg 1998; Tsuchiya 2002), 
and the third is an impulsive deposition of momentum to the gas disk by passage through it 
(Ibata & Razoumov 1998). The direct tidal torque for a satellite of mass m and distance r 
scales as m/r 3 . Therefore, the direct tidal effect of Sgr is no stronger than that of the LMC, 
whose direct tidal torque is not sufficient to induce the warp (Hunter & Toomre 1969). This 
means that the gravitational torque of Sgr itself cannot be the coupling mechanism. 

If the primary perturber is instead a wake in the halo, the strength of the torque scales 
as ^wakc/ r wakc- The mass of the wake scales as the mass of the satellite and as the density 
of the halo at the wake radius (Weinberg 1998). The wake develops at half the satellite's 
orbital radius (Tsuchiya 2002). Therefore, for an isothermal halo, the strength of the torque 
scales as m/r 5 . In this case, the effect of Sagittarius is 10-50 times stronger than that of the 
LMC. It is plausible that in Tsuchiya (2002) 's lower mass simulation, in which the LMC did 
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not excite a warp, a satellite with Sagittarius's parameters would have. Further simulations 
which better reproduce the observed Sgr-M W system could confirm or falsify this suggestion. 

Ibata & Razoumov (1998) suggest that the impulsive deposition of momentum to the 
gas disk could excite the warp. The mass they use for Sgr, 5 x 10 9 M , is quite large, and they 
find very little warping in their 1 x 10 9 M simulation. However, they only model a single 
interaction. In order for the angular momenta to reach an equilibrium, as they appear to have 
done, there must be repeated or continual encounters. Helmi & White (2001) find orbital 
periods of ~ 1 Gyr for Sagittarius, indicating that it has passed through the disk several 
times. Further simulations that follow the evolution of the system over many encounters are 
necessary to better understand the predictions of this model; meanwhile, it cannot be ruled 
out. 



5. Summary 

The orbital angular momentum of the Sagittarius dwarf galaxy and the component of 
the Milky Way disk angular momentum due to the Galactic warp are both directed toward 
/ ps 270°, 6 = 0° with magnitude 2-8 x 10 12 M Q kpc km s _1 . Such a coincidence suggests that 
they are a coupled system, i.e., that Sgr is responsible for the warp. The direct gravitational 
tidal torque of Sgr cannot cause the warp. Interaction via a gravitational wake in the Milky 
Way's dark matter halo, and impulsive deposition of momentum into the disk by passing 
through it are still both possible coupling mechanisms. More simulations of each of these 
models are necessary to discriminate between their effects. 

Many thanks to Casey Meakin for useful discussions and comments. 
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